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Abstract. 

We present a simple and thermodynamically consistent cosmology with a 
phcnomenological model of quantum creation of radiation due to vacuum decay. 
Thermodynamics and Einstein's equations lead to an equation in which H is 
determined by the particle number N . The model is completed by specifying the 
particle creation rate T — N/N, which leads to a second-order evolution equation 
for H. We propose a simple T that is naturally defined and that conforms to the 
thcrmodynamical conditions: (a) the entropy production rate starts at a maximum; 
(b) the initial vacuum (for radiation) is a non-singular regular vacuum; and (c) the 
creation rate is initially higher than the expansion rate H, but then falls below H . 
The evolution equation for H then has a remarkably simple exact solution, in which a 
non-adiabatic inflationary era exits smoothly to the radiation era, without a reheating 
transition. For this solution, we give exact expressions for the cosmic scale factor, 
energy density of radiation and vacuum, temperature, entropy and super- horizon scalar 
perturbations. 



PACS numbers: 9880C, 9880H, 9530T, 0540 



1. Introduction 

The creation of matter and entropy from vacuum has been studied via quantum field 
theory in curved spacetime (see for example 0]). Most cosmo logical models exhibit 
a singularity which presents difficulties for interpreting quantum effects, because all 
macroscopic parameters of created particles are infinite there. This leads to the problem 
of the initial vacuum. A regular vacuum for a species of created particles can be defined 
in simple terms as a state where all mean values describing the particles, such as energy 
density, number density, entropy etc., are zero. But this simple condition is not achieved 
in many scenarios, so that either one has to postulate an initial state beyond the 
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singularity, or to assume that there was a nonzero number of particles at the initial 
vacuum. 

One attempt to overcome these problems is via incorporating the effect of particle 
creation into Einstein's field equations. For example, in the papers of the Brussels 
group [[J, the quantum effect of particle creation is considered in the context of the 
thermodynamics of open systems, where it is interpreted as an additional negative 
pressure, which emerges from a re-interpretation of the energy-momentum tensor. This 
effect is irreversible in the sense that spacetime can produce matter, leading to growth 
of entropy, while the reverse process is thermodynamically forbidden. 

The main difference with our present paper is that in || the law of massive particle 
creation, i.e. the mechanism of energy flow from the gravitational field to matter, 
leads to a non-zero number of particles at the beginning of expansion, described as a 
fluctuation of the regular vacuum. These results were recently generalized in a covariant 
form in j|. Our approach differs from that of 0, Q in that we do not modify the field 
equations. Instead, we interpret the source of created particles as a decaying vacuum, 
described phenomenologically by a time-dependent cosmological 'constant' A(t). 

A number of decaying vacuum models has appeared in the literature (see 0, |J 
and references cited there). Inflationary models with fixed cosmological constant and 
cold dark matter have been successful in accounting for the microwave background and 
large-scale structure observations, while also solving the age problem. However, these 
models are challenged by the reduced upper limits on A arising from the Supernova 
Cosmology Project, and also by the long-standing problem of reconciling the very large 
early-universe vacuum energy density with the very low late-universe limits ||]. One 
resolution of these problems is a decaying A. In common with ||, we attempt to 
provide some clear and consistent physical motivation for the particular form of vacuum 
decay, rather than an ad hoc prescription. In ad hoc prescriptions, the functional form 
of A(t) or A (a) (where a is the scale factor) is effectively assumed a priori. Often 
power-law forms for A are assumed (see, for example, and references cited there). 
Exponential decay laws have also been assumed || . Typically, the solutions arising from 
ad hoc prescriptions for A are rather complicated, and moreover, it is often difficult to 
provide a consistent simple interpretation of the features of particle creation, entropy 
and thermodynamics. 

In contrast to many other models, we propose a simple, exact and thermodynami- 
cally consistent cosmological history. The latter originates from a regular initial vacuum 
with a maximal initial entropy production rate. Together with a naturally defined cre- 
ation rate, this leads to a simple expansion law and thermodynamic properties, and to 
a definite estimate for the total entropy in the universe. The very existence of an initial 
maximal entropy production rate reflects a subtle interplay between the conservation 
equation, the second law of thermodynamics and Einstein's equations, which is at the 
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heart of our model. 

Non-adiabatic inflationary models differ from the standard models (see for example 
0), in that: (a) radiation is created continuously during inflation, rather than during 
reheating; (b) the continuous vacuum decay itself initiates a smooth exit from inflation 
to the radiation era; (c) entropy and heat production take place continuously, without 
the need for reheating. In the standard approach, the scalar field drives adiabatic (i.e., 
isentropic) inflation, followed by a non-equilibrium reheating era when the field decays 
into radiation and inflation is brought to an end. The potential of the field is then the 
key ingredient. In the alternative approach, the key ingredient is essentially the model 
of vacuum decay. In contrast to the ad hoc models that assume functional forms for A, 
and as an alternative to the thermodynamically motivated model of f|, we arrive at a 
phenomenological model of decay by imposing simple and thermodynamically consistent 
physical requirements. A related physically consistent model can be found in the 'warm' 
inflationary scenario ]lJJ . 

The first law of thermodynamics for open systems (with particle creation) and 
Einstein's equations lead to a first-order equation for the expansion rate H = a/ a, 
whose source term is determined by the particle number N. A further equation in H 
and N arises from a simple model for the particle creation rate T = N/N. We impose the 
thermodynamical non-equilibrium condition that T, and hence the entropy production 
rate, is maximal at the beginning of expansion. A further initial condition is that the 
initial vacuum for radiation is non-singular. We also require that T > H initially, so that 
created particles are thermalized, while T < H later on, as particle production becomes 
insignificant. Our simple model for T is naturally defined by the gravitational dynamics 
and conforms to these thermodynamical requirements. We decouple the equations to 
get a second-order evolution for H, and we find a remarkably simple exact solution 
H(a). 

Since the exit from inflation to the radiation era is smooth, we avoid the problem 



of matching at the transition. A similar smooth evolution has been used in [IT], [12], [TJ 
but in the context of adiabatic inflation, and without a consistent physical foundation. 



In effect, we show that the ad hoc form of H(a) given in |TTJ follows from our simple 



physical conditions and thermodynamic arguments. In [14], a kinematic analysis is given 
for various non-adiabatic inflationary evolutions with smooth exit, but these evolutions 
are outside the scope of our model. 

The choice of a as dynamical variable and the very simple form of H (a) that meets 
the physical conditions, lead to elegant expressions for all parameters describing the 
radiation and decaying vacuum, and also to a physically transparent interpretation of 
these results, including the estimate of entropy. In addition, the equation for super- 
horizon scalar perturbations can be solved exactly for this form of H(a). 

We present in Sec. 2 the evolution equation for H(a) and its simple solution, that 
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follow from our simple physical constraints. In Sec. 3 we analyze the thermodynamics 
of the radiation produced in the course of vacuum decay, and estimate the entropy of the 
created radiation. Sec. 4 contains a summary and concluding remarks. In a subsequent 
paper, we will discuss generalizations of the present model. 
We use units with StcG, c and k& equal 1. 



2. The simple model 

Consider a spatially flat Friedmann-Lemaitre-Robertson- Walker universe 

ds 2 = g^dx^dx" = dt 2 - a 2 (t) dx 2 + dy 2 + dz 2 , 
containing matter with equation of state 

P = 7P, 

where 7 (0 < 7 < 1) is a constant parameter. Later we will specialize to the case of 
radiation, 7 = |. The energy-momentum tensor of matter is 

T fu = p(t) [(7 + i)w M Wi/ + ig^] , = 5^ , 

while the energy-momentum tensor corresponding to the quantum vacuum energy is 

rfi = ($$,) = A(t)g, u . 
Then the conservation equations V V {T^ + T®) = reduce to 

p + 3( 7 + l)tfp = -A, (1) 

showing how energy is transferred from the decaying vacuum to matter. Note that (|l|) 
is equivalent to the energy balance of an imperfect fluid with scalar viscous pressure 

This is an example of the known result that cosmological particle production may be 
interpreted as an effective bulk viscous pressure (see [15) and references cited there). 
The field equations R^ u — \Rg tlv = Tj^ + T^ u are 

3H 2 =p + A, (2) 

2H + 3H 2 = - 7p + A , (3) 

and if both are satisfied then the energy equation ([I]) follows identically. 

Following [11], we use a as a dynamical variable instead of t, and so we consider the 
Hubble rate as H = H(a). Then equations (2) and (3) give 

p(a) = - aH(a)H>(a) , (4) 

A(o) = 3H 2 {a) + (—-fj aH(a)H'{a) , (5) 
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where primes denote d/da. Given H(a) and 7, we can calculate p(a) and A (a). In order 
to determine further properties of H(a), we will impose simple physical requirements. 

We assume that the evolution of p is governed not only by expansion but also by 
the creation of particles, i.e. from the thermodynamic point of view, we have an open 
system. According to || , the first law of thermodynamics generalized for open systems 
is 

d(pV) + pdV - ( d(nV) = , (6) 



n 

where p = 7/? is the pressure, n = N/V is the particle number density, A" is the number 
of particles in the observable universe, and V oc a 3 is the comoving volume of the 
observable universe. The thermodynamic equation (^) implies 

(pa 3 )' (a 3 )' ,,,.N' 
which integrates to 

A (NV I+1 

9 = (7TT) taO ' (7) 

where A is a positive constant. From now on, we take 7 = |, i.e. we assume that only 
radiation is created. This will also allow us to define entropy production via the photon 
number. 

By equation (4), equation (0) can be rewritten as an evolution equation for H , with 
source term determined by the particle number: 

This is the fundamental equation which follows from Einstein's equations and the 
thermodynamic equation (|5]). The creation mechanism is phenomenologically encoded 
in the source term N(a). The model must be completed via an equation that determines 
N. If the creation rate of radiation r = N/N = aHN'/N is given, then equation (§) 
implies the second-order equation 

3aHH" + 3aH' 2 + (15H - AT) H' = 0. (9) 

We seek a model in which most of the particle creation effectively takes place in the very 
early universe, starting from a regular vacuum. More precisely, we impose the following 
thermodynamical requirements: 

(a) Maximal entropy production rate (equivalently, maximal particle creation rate) at 
the beginning of expansion, so that the universe starts in a state furthest away from 
equilibrium and then tends toward equilibrium as the expansion proceeds. 

(b) A true (regular) vacuum for radiation initially, so that p — ► as a — > 0. 

(c) T > H in the very early universe, so that we can treat the created radiation as 
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forming a thermalized heat bath. Subsequently, the creation rate should fall behind the 
expansion rate as particle creation becomes dynamically insignificant. 

The fundamental physical quantities that are naturally defined by the gravitational 
dynamics in our model are the expansion rate H and the total energy density U = p + A. 
Both of these quantities can define in a natural way a gravitational creation rate T. The 
simplest model T oc H fails to satisfy requirement (c) above. Furthermore, it leads to a 
solution H(a) of the evolution equation (Q) which violates requirement (b). Therefore 
we propose r oc U. This will satisfy requirement (a) if (b) holds, since the initial 
condition (b) implies that U approaches its maximum value U(0) = A(0) as a —>■ 0. 
Below we show that the solution of equation (§) does verify condition (b). The increase 
of p due to creation in the very early universe partially offsets the decrease in A, so 
that U decreases slowly in the very early universe, and the entropy production rate 
remains high. Requirement (c) will be satisfied because in the radiation era, A becomes 
negligible and p decays like H 2 . Thus T decreases more rapidly than H, and will have 
fallen below H at some epoch. 

Hence we propose the simple model that the particle creation rate is proportional to 
the total energy density. By the Friedmann equation (2), it follows that 

r=aHc {§~y !) (io) 

where a is a dimensionless free parameter, and H e = H(a e ), where a e is some fixed 
epoch. Equation (0) becomes the decoupled, second-order evolution equation for H (a) : 

3aHH" + 3aH' 2 + 15HH' - A—H 2 H' = . (11) 



This equation has the first integral 



Aq, 

3aHH' + 6H 2 - — -H 3 = , (12) 
3H 



1 c 



where we have used the fact that H and aHH' tend to zero for large a (i.e. in the 
radiation era), in order to remove a constant of integration. The solution of equation 

(0) ^ 

1/2 

3' 



_9_ /He 

2a \H 



a 

Me 

where /3 is a constant, and we have taken into account that if is a decreasing function. 
Evaluating at a = a e , we see that 

l + ,3 2 = ^, (13) 

which implies the constraint a < | on the creation parameter a. Re- arranging the 
solution, we find the remarkably simple form for the expansion rate that follows from 
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our thermodynamic model: 



f = (l + /5 2 ) 



a e 2 + /3 2 a 2 



(14) 



This solution approaches de Sitter inflation as a — > 0, i.e. H ~ constant, and it 
becomes radiation-like for a — > oo, i.e. if ~ a~ 2 . It follows that as a —>■ 0, the cosmic 
proper time t — > — oo. A naturally defined epoch is a ex = a{t ex ) of exit from inflation, 
which is defined by a(t ex ) = 0, or equivalently H(a ex ) = —a ex H'(a ex ). It follows from 
(|H|) that a e = fta ex . 

These results reflect a subtle interplay between Einstein's gravitational dynamics 
and thermodynamical constraints on particle production. The decay of the vacuum 
into radiation drives inflation, but the same decay, by reducing the vacuum energy 
density, leads to a deceleration of expansion and a smooth exit from inflation. Since 
the initial vacuum for radiation is regular, the particle number and hence entropy are 
initially zero. As we show below, the total entropy produced in the observable universe 
in an infinite time is finite. This is consistent with the existence of a smooth exit, since 
unending inflation would produce infinite entropy. The initial rate of entropy production 
is maximal, reflecting the feature that the universe starts furthest from equilibrium 
and approaches asymptotically a state of equilibrium. We note also that the initial 
entropy production rate is a finite maximum value, rather than being unbounded from 
above. The latter possibility is ruled out by the de Sitter-like nature of the inflationary 
expansion, which implies that the expansion rate H is bounded from above (unlike 
power-law inflation for example). 

The freedom in ft, or equivalently a, by equation (|T3|), provides us with an extra 
adjustable parameter. However, for simplicity, we will not use this freedom, since the 
subsequent results are not modified in any essential way for general ft. Henceforth we 
take ft = 1, i.e. a = |, which means that a e = a ex . Thus, we arrive finally at the simple 
expansion rate 



H(a) = 2H C 



a 2 + a 2 



(15) 



This form of H(a) was given in [TT| as an ad hoc toy model to achieve smooth exit 
from inflation to radiation, but without a physical basis such as that given here. The 
expression for the cosmic proper time follows on integrating equation flTB]): 



t = t e + 



In ( — 

a e 



a \ 2 / a x 2 



-) -1 



3. Thermodynamics of radiation 

On substituting equation flTSP into equations (4) and (5), we find exact expressions for 



the energy density of radiation and the vacuum: 



p(a) = 12HI ( f^y 3 , A(a) = 12H> f-^V • (16) 



It follows that A(0) = 12F e 2 . Note that (|T|) implies H(0) = 2H C . Note also that the 
effective bulk viscous pressure arising from particle production has the form 
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3H \9HJ y \a 2 p + a 



Now p reaches a maximum at a m = a c /V% with 

Pm = p(ctm) =f ff e 2 , A(a m ) = 2p m . 
Note also that p and A are equal at exit: 

p(a e ) = A(Oe) = f^e > 
while for a ^> a c , i.e. during radiation-domination, 

p(a) ~ — , A(a) ~ — , 
a 4 a 

so that A rapidly becomes negligible. 

The formulas (|16|) reflect the creation of radiation due to vacuum decay. The initial 

value p(0) = confirms that the field corresponding to radiation is initially in a regular 

vacuum state. This means that the formula ( |TE| ) gives an absolute measure of radiation 

produced in the universe, as a result of conversion of energy from the vacuum described 

by A. 

Substituting equation ( |15|) into equation (^) we get the exact form for the particle 



number 

2 \ 9/ 4 

al + a 2 



N{a) = N 00 [^—\ , (17) 



where iVoo is usually taken to be about 10 88 . It follows that about 2 x 10 87 particles 
have been created at exit. Note that the initial conditions and the evolution equations in 
our model imply that a finite number of particles is produced in the observable universe 
during the entire expansion. 

Since N(0) = n(0) = 0, the initial state of the field has no particles, i.e. it is a 
regular vacuum. The number density is 

"W^M^H^f- (18) 
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and reaches its maximum also at a m . 

The creation rate of radiation is given by equations (|10|) and (|T5|) as 

' " 2 N * 9,. /H\ 2 



r = 9ft^j =- 4 HhJ ■ < 19) 

In order to define the radiation temperature (and then entropy), we would like to invoke 
the standard black-body relation. A justification of this is as follows. From equation 
([19]), it follows that the creation rate T(a) exceeds the expansion rate if (a) for 

a < J\ a e . 

Thus it is reasonable to treat the created radiation as forming a thermalized heat bath 
in the initial stage of expansion. For a > ^J\a e , the creation rate falls behind the 
expansion rate, and created particles will be out of equilibrium. However, by this stage 
of the expansion, the energy density in newly created particles is too small to disturb 
the effective thermalization. Thus it seems reasonable to use the black-body relation 
for the radiation throughout the expansion, and to define the temperature by 

3n(a) 2 i / 4 n e \ aj \a£ + a 1 J 

where we have used equations (|T6|) and (|T8|) . At the initial radiation vacuum, it is clear 
that T(0) = 0, and T increases to its maximum value at a m : 

.27/ n P ' 



T m = T(a rj 



This temperature may be thought of as analogous to the reheating temperature in 
standard models. During the radiation era, i.e. for a ^> a c , 

T ~ a -1 , 

in agreement with the standard result for free radiation in an expanding universe. 
The formulas for p and n can be presented in the thermodynamic form 

, = 24(1)7-, - 8 (S) r3 ' P1) 
Combining now the Gibbs equation 
TdS = d(pV) + pdV , 

with equation (|6]), and using the definition pop of T, we obtain the entropy of radiation 
in the observable universe as 



S(a) = 4N{a) , 



(22) 
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so that 5(0) = as expected. This gives a reasonable value for the entropy produced 
during the overall evolution of the universe: 

5oo = 4^0,, « 4 x 10 88 . (23) 

In the standard model (adiabatic inflation followed by reheating), one can estimate the 
entropy production due to reheating by matching exact de Sitter inflation to an exact 



radiation era, with an instantaneous transition at a e . This gives |Ll 

C _ 4 1/4 3/4 _ o 
°e — 3 y He — °oo i 

where g ~ 100, and leads to a value of the same order of magnitude as our result. 
4. Conclusion 

We have considered a simple and thermodynamically consistent scenario encompassing 
the decay of the vacuum, the creation of radiation and entropy, and a natural smooth 
transition from inflationary to radiation-dominated expansion. In order to treat all 
matter in the universe as created from a regular vacuum, we impose the condition 
that p — > as a — > 0. We impose the further initial condition that the entropy 
production rate is a maximum. A simple model for the particle creation rate T, i.e. 
that T is proportional to the comoving total energy density, is shown to be consistent 
with these initial conditions and with the requirement T should start above and then 
fall below the expansion rate H. We showed that the field equations and the first law of 
thermodynamics @ generalized for open systems with creation of matter, then imply 
a second-order evolution equation ( |Tl~D for H(a). This equation has the remarkably 
simple solution (|i~4), and we used this together with black-body thermodynamics and 



the Gibbs equation to define the temperature and entropy. 

Our postulate T oc H 2 for the creation rate, can be contrasted with other work. 
For example, in ||, N oc H 2 , while T is constant; in 0, A oc exp(— t/r), and there 
is no simple expression for T in terms of H. The postulates in these papers are ad 
hoc, whereas we have tried to give a thermodynamic justification for our postulate. It 
is also interesting to compare our decaying- vacuum/ radiation model with scalar-field/ 
radiation models |T^, |TH]. In these latter models, a phenomenological term is introduced 



into the Klein-Gordon equation for the scalar field 0, describing the interaction between 
the decaying field and radiation. The resulting energy balance equation is 

p + 4Hp = T^ 2 , 

where is the phenomenological decay rate. Comparing this with our energy balance 
(|l|), and using (0), we see that T^cf 1 corresponds to our |Tp. It is not obvious whether 
a scalar field interpretation of our model exists that would produce the rate oc H 2 . 
This is a subject for further investigation. 
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An important feature of the Hubble rate flT5|) is that lim a ^ H(a) = constant and 
lim a ^ H'(a) = 0. According to equations ([K]) and ([T7|), this avoids any divergences in p 
and N at a = 0. In addition, it follows from the equations (|16D, (|TT1), fl20|) and (^) that 
all thermodynamic parameters describing the created radiation vanish initially, i.e. all 
forms of matter which we observe now in the universe were in a state of regular vacuum, 
and all constituents of the universe were created from this vacuum in the course of the 
decay of the vacuum energy density A from its initial value 12H^. We also showed 
that our requirements imply finite particle and entropy production during the entire 
expansion of the observable universe. 

Our simple model leads to exact expressions for t(a) and for the thermodynamic 
variables. In addition, the form ( |i~5l) of H(a) leads to an exact solution for the large- 
scale modes of scalar perturbations, described gauge-invariantly by Bardeen's potential 
$. The solution is |17|: 

where C' ± = 0. We will not discuss the important question of the source of these 
perturbations, except to mention that, as pointed out in |I(J, non-adiabatic inflationary 
scenarios allow for the possibility that the seeds of density perturbations are of 
predominantly thermal, rather than quantum, origin. 

Further discussion of the physical processes taking place during non-adiabatic 
inflation, and of their effects on baryogenesis, the microwave background radiation and 
structure formation, is given in |10|, [14], |7], |], [18| . 
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